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ABSTRACT 

We present a detailed analysis of 17, 852 quiescent, Luminous Red Galaxies 
(LRGs) selected from Sloan Digital Sky Survey (SDSS) Data Release Seven (DR7) 
spanning a redshift range of 0.0 < z < 0.4. These galaxies are co-added into four equal 
bins of velocity dispersion and luminosity to produce high signal-to-noise spectra 
(> 100A -1 ), thus facilitating accurate measurements of the standard Lick absorption- 
line indices. In particular, we have carefully corrected and calibrated these indices onto 
the commonly used Lick/IDS system, thus allowing us to compare these data with 
other measurements in the literature, and derive realistic ages, metallicities ([Z/H]) 
and a-element abundance ratios ([a/Fe]) for these galaxies using Simple Stellar Pop- 
ulation (SSP) models. We use these data to study the relationship of these galaxy 
parameters with redshift, and find little evidence for evolution in metallicity or a— 
elements (especially for our intermediate mass samples). This demonstrates that our 
subsamples are consistent with pure passive evolving (i.e. no chemical evolution) and 
represent a homogeneous population over this redshift range. We also present the 
age-redshift relation for these LRGs and clearly see a decrease in their age with red- 
shift (~ 5Gyrs over the redshift range studied here) which is fully consistent with the 
cosmological lookback times in a concordance ACDM universe. We also see that our 
most massive sample of LRGs is the youngest compared to the lower mass galaxies. 
We provide these data now to help future cosmological and galaxy evolution studies 
of LRGs, and provide in the appendices of this paper the required methodology and 
information to calibrate SDSS spectra onto the Lick/IDS system. 

Key words: methods: observational - galaxies: elliptical - galaxies: abundances - 
galaxies: evolution - galaxies: fundamental parameters 



1 INTRODUCTION 

The formation and evolution of massive elliptical (passive) 
galaxies in the Universe is interesting for both studies of 
galaxy evolution and cosmology. In the former case, such 
galaxies present an observational challenge for hierarchi- 
cal models of structure formation, as some form of feed- 
back is required to suppre ss on-going star f ormation in 
such massive systems (see IWake et al.l (|2006l ) and refer- 
ences therein). For cosmology, such massive ellipticals can 
be used to directly constrain cosmological parameters (e.g., 
as standard candles: ICollins fc Mannl (| 19981 )') and provide ef- 
ficient tra£e£s_of_th*I^r^eScale Structure (LSS) in the Uni- 
verse (fEiscnste in et al]|200ll ). Moreover, the (relative) ages 
of massive ellipticals, as a function of redshift, offers the pos- 
sibility of directly constraining the Hubble parameter, thus 
providing vital information on the expansion history of the 
Universe and therefore, the equation of state of dark en- 



ergy (see ? for a disc u ssion of the u nderly ing c oncept, and 
IJimenez et~ai] (120031 ): ISimon et all (|2005h and IStern et al] 
|2009l) for observational constraints obtained from using this 
technique) . 

In this paper, we revisit the techniques used to constrain 
cosmological parameters through the age-redshift relation- 
ship of passive (elliptical), massive galaxies. In detail, we 
present a n ew analysis of the ages of Luminous Red Galax- 
ies (LRGs; lEisenstein et alj (|2001)), as selec t ed fro m the 
Sloan Digital Sky Survey (SDSS: lYork et all (|2000h ). The 
key difference present herein co mpared to other such w ork 
with SDSS spectral d a ta (e.g. IBernardi et ail (|2003al lbh: 
lEisenstein et al.l <|2003h ; IBernardi et al l i|2006)) is that we 
calibrate, for the first time, the SDSS spectr a onto the well- 
known , and well-studied, Lic k/IDS system (jBurstein et alj 
(QUI); IWorthev et~ail l| 19941 )), thus allowing us to exploit 
a host of previous work on this system including the Sim- 
ple Stellar Population (SSP) modeling in the literature 



Dan P. Carson and Robert C. Nichol 



(|Worthevl (|l998l 'l; lMarastonl (|l998l l: iKorn et all |2005l )l and 
comparisons with Lick measurements in the literature. In 
a companion paper, we will use the age-redshift data de- 
rived for SDSS LRGs in this paper to obtain cosmological 
constraints. 

In Section 2 of this paper, we outline the SDSS galaxy 
data we use to construct our age-redshift relation, specifi- 
cally luminous (massive), quiescent galaxies. We also discuss 
the stellar data available to calibrate these galaxies onto the 
Lick/IDS system (see also Appendix A). In Section 3, we 
provide details about the SSP models we have used to de- 
termine the ages, metallicities and a-enhancements of our 
galaxies, using a Markov Chain Monte Carlo (MCMC) tech- 
nique for parameter estimation. We discuss our results in 
Section 4, including the effect of priors, and conclude in Sec- 
tion 5. Appendix A provides the extensive details on how 
we calibrate the SDSS spectra onto the Lick/IDS system, 
including all the necessary corrections made to the data. In 
Appendix B, we provide a review of tests we have performed 
to quantify the robustness of our spectral measurements, in- 
cluding comparisons with data from the literature. We as- 
sume a concordance, flat ACDM cosmology, with h = 0.7 
and A = 0.7, where required. 



2 SAMPLE SELECTION 



The Sloan Digital Sky Survey (SPSS; lYork et~all |2000| ; 
IStoughton et al.ll2002l ; lAbazaiian et aLlbOol T is a photomet- 
ric and spectroscopic survey, covering ~ 8000deg 2 of the 
northern sky, using a 2.5 meter telescope at Apache Point 
Observatory in New Mexico, USA. The photometric survey 
consists of simultaneous observations of the sky using 5 opti- 
cal filters (u, g, r, i and z), providing a database of hundreds 
of millions of detected object s all with accurate photometric 
and astrometric calibra tions (|Pier et al 1 l2003l ; lLupton et all 
l200ll ; lHogg et al.ll200ll l. 

The SDSS spectroscopic galaxy survey consists of two 
samples of galaxie s selected using diff erent criteria; namely 
the MAIN sample jStrauss et al.ll2002T ) and the LRG sample 
l|Eisenstein et al.ll2001 ). The SDSS spectra of these galaxies 
span a wavelength range of 3800 < A < 9200A with a median 
resolution of R ~ 1800 (R = A/ A A), which is approximately 
3 A (although this varies as a function of wavelength and is 
different for each fiber). The SDSS spectra are automatically 
reduced using dedicated software, which flux calibrates the 
spectra and references them to the heliocentric frame and 
converts to vacuum wavelengths. The software also mea- 
sures a redshift for each object and measures a series of 
spectral features consistent (in their wavelength definitions) 
as the standard Lick indices (see Section 2.2), but are not 
calibrated onto the Lick systems as no attempt has been 
made to match to the resolution of the Lick/IDS system. 

For the analysis presented in this paper, we do not 
use the Lick indices measurements from the standard SDSS 
pipeline, but instead determine our own line-strengths after 
matching the instrumental resolutions (see Section[A]). How- 
ever, we do use redshifts (z), velocity dispersions (velDisp), 
magnitudes (u,g,r,i,z) and other derived quantities such 
as the r-band de Vaucouleurs radii (deVRacLr) and the r- 
band de Vaucouleurs profile axis ratio (deVAELr) from the 
standard SDSS spectral pipeline, made available through 



the Catalog Archive Server (CASfl The spectral data used 
in this paper was obtained from the Data Archive Server 
(DASfl 



2.1 SDSS Luminous Red Galaxies 

For our analysis, we only used Cut I Luminous Re d 
Galaxies (LRGs) as outlined in (|Eisenstein et al.1 1200 ll ). 
This is achieved by selecting objects from CAS using the 
TARGET_GALAXY_RED flag, thus yielding a pseudo volume- 
limited sample of LRGs between 0.15 < z < 0.35. Below 
z = 0.15 contamination by low redshift star-forming galaxies 
increases the space density of galaxies that satisfy the SDSS 
LRG colour-colour selection, while above z = 0.35 the space 
density of SDSS LRGs decreases due to the 4000 A break 
dropping out of the SDSS gr-band. We do not use the Cut 
II LRG sample at z > 0.4 as we require velocity dispersion 
measurements for each object and these are not supplied 
above this redshift. Therefore, we impose an upper redshift 
limit of z — 0.4 in our sample and analysis. 

We further restrict our sample to only include LRGs 
with specClass EQ ' SPEC_GALAXY' (spectrum classified as 
a galaxy), zStat EQ 'XC0RR_HIC (redshift obtained from 
cross-correlation with a template), zWarning EQ (no 
warning flags on the redshift), eClass < (indicates an 
old stellar population), z < 0.4 (redshift less than 0.4) 
and f racDev_r > 0.8 (indicating a surface brightness profile 
best fit by a de Vaucouleurs profile). Based upon these se- 
lection criteria, we obtain a sample of approximately 77, 000 
LRGs (see Figure [2] for the redshift distribution of this sam- 
pie). 



2.1.1 Quiescent Galaxies 

At low redshift, our sample will be contaminated by bulges 
in late-type galaxies due to the fiber size (3") of the SDSS 
spectrographs. To reduce this contamination, and increase 
the number of truly quiescent galaxies in our sample, we 
make further cuts based on the emission line properties of 
our LRG sample, e.g., we u se standard emission l ines such as 
Ha, Hff and [OHI1A5007 dFukugita et all |2004| ; [ Zhao et al.l 
2006; Rose boom et al. 2006). If left unchecked, such nebu- 
lar emission, or emission associated with low-ionization nu- 
clear emission-line region (LINER) activity, would confuse 
our interpretation of the SSP parameters derived from these 
objects. 

To combat this, we use the MPA-JHU spectral line 
dataS to select only those objects that are consistent 
with z ero emission. In their original work, iTremonti et al.l 
||2004|) fitted the sp e ctral energy distributions (SEDs) 
of iBruzual fc Chariot! (|2003l ) to the SDSS galaxy spec- 
tra in Data Release Four (DR4), subtracted off the best 
fitting SED, and then measured the line-strengths of 
[OII]A alr 3726, H/3 , [O III]A ai r5007, [N II]A Qlr 6584, Ha and 
[S II]A air 6717. This SED fitting approach has now been ap- 
plied to the SDSS Data Release Seven (DR7) dataset using 



1 http:/ /www.sdss.org/dr7/access/index.html#:CAS 

2 http:/ /www.sdss.org/dr7/access/index.html#:DAS 

3 http:/ /www. mpa-garching.mpg.de/SDSS/DR7/raw_data. html 
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Figure 1. Distribution of the line-strengths from the MPA-JHU dataset for the LRGs of our sample. Dotted contours show the full LRG 
sample while solid lines show the quiescent sample selected based on their Hce and [Oil] line-strengths. While only two lines have been 
used in the production of the quiescent sample, all lines show distributions consistent with zero-emission. 



the stellar population synthesis spectra of Chariot & Bruzual 
(2008). We use this latest DR7 dataset her e in. 

In an approach similar to lYan et alJ (|2006T) we fit a 
t wo component fu nction to the line-strength distributions 
of iTremonti et all that consists of a normal component to 
describe the quiescent LRGs, plus a log-normal component 
to describe the active objects. The best fit is then used to de- 
te rmine the small zero-point offsets that exist in the dataset 
of ITremonti et al.l from which we then select objects that 
are consistent with zero emission, at the 2cr level, in Ha and 
[O II] , hence removing the need to make corrections to the 
Lick index line-strengths for the possible presence of neb- 
ular emission. With these constraints on the emission-line 
characteristics of the LRG sample, we obtain approximately 
35, 000 galaxies that we refer to as the quiescent LRG sam- 
ple. 

While it is possible to use different, and more, emission 
lines in determining the quiescent sample, doing so results 
in a significant reduction of our sample size. For instance, 
including H/3and [O III] in the joint constraint with Ha and 
[O II] yields approximately 28, 000 LRGs over 0.0 < z < 0.4 
and further inclusion of [N II] and [S II] constraints yields 
approximately 19, 000 objects. While the more restrictive 
sample selection ensures a more quiescent sample, the re- 
duction in the total number of objects available hampers 
our ability to explore the evolution of the SSP parameters 
over the desired redshift range. But, even selecting on Ha 
and [O II] alone helps to ensure that we are not including 
objects with significant emission in H/3 as show in Figure [T] 
In Figure we also show the redshift distribution of the 
quiescent LRGs sample. 



2.1.2 Further Selection Criteria 

In order to select the same population of objects over 
the redshift range 0.0 < z < 0.4 we need to select 
on physical properties such as velocity dispersion and 
absolute luminosity. Aperture corrections are applied to 
the SDSS velo city dispersion measur ement s following the 
approa ches of |j0rgensen et al.l (| 19951 ^) and iBernardi et al.l 
(2003a), see Section IA5I To select on luminosity we have 
performed K-corrections using the code kcorrect code of 
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Figure 2. Redshift distribution of our LRG sample (solid line) 
and our quiescent sample (dashed line). 

Blanton & Roweisl (|2007f ) and have K-corrected all galax- 
ies to a redshift of 0.2 (the median of our distribution) in 
order to minimise errors in these K-corrections. 

We have modeled the luminosity evolution in our LRGs 
using a simple linear fit to the volume limited part of the K- 
corrected, magnitude-redshift distribution and then extrap- 
olated this to higher and lower redshifts, see Figure [3] From 
this, we obtain an evolutionary (e) correction which we ap- 
ply to our quiescent LRG sample. A more detailed approach 
to determine the evolution correction using stellar popula- 
tion synthesis models is unlikely to significantly improve the 
accuracy of the e-correction given the difficulty that current 
models have in explaining the evolution of the S PSS LRGs 
colours (|Wake et alj|200a iMaraston et al.ll2009T ). 

When reporting absolute magnitudes, we will use the 
notation Mo.2 r to denote r-band SDSS magnitudes that have 
been K+e corrected such that the r-band has been shifted 
to z = 0.2. 

After correcting velocity dispersions for aperture effects 
and performing K+e corrections to the magnitudes, we have 
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Figure 3. Distribution of our quiescent LRG sample in the 
redshift-luminosity plane, with a simple linear model (green/solid 
line) used to describe the luminosity evolution of the LRGs. Con- 
tours show the distribution in number density of the whole sam- 
ple, and are in the intervals 5,10,50,100. The region bounded by 
the vertical dotted lines represent the volume-limited portion of 
our sample and the (red/dashed line) shows the median luminos- 
ity in redshift intervals of Az = 0.01. 



produced four subsamples of quiescent LRGs (see Table 2 for 
numbers in each sample). These subsamples span an inter- 
vals of AMo.2 r = 0.6 and Aa v = 30kms _1 and are shown as 
the boxed regions in FigureU] Objects satisfying these crite- 
ria are then co-added to produce high signal-to-noise spectra 
after binning into redshift intervals of Az — 0.02. This red- 
shift interval is small enough that there is little evolution 
in the ages of the objects within the bin such that we are 
co-adding objects of sufficiently similar ages, e.g. in the case 
of a concordance cosmology the age evolution within a bin is 
approximately 0.2 Gyr over the redshift range 0.0 < z < 0.4, 
and is big enough to yield a sufficient number of objects to 
reach S/N levels of approximately 100 A -1 

In generating the co-added spectrum we have adopted 
an approach similar to lBernardi et all (|2006l ). Briefly, after 
binning our sample in redshift, absolute magnitude and ve- 
locity dispersion we shift each spectrum to its restframe. All 
spectra in a given bin are normalised to the median flux in 
the region 4500 — 5500 A and combined pixel-by-pixel using 
a weighted arithmetic mean where the weights are deter- 
mined by the inverse variance of the flux in each pixel. This 
approach means that not all spectra contribute equally to 
the co-added spectrum but it does mean that those pixels 
contaminated by residuals from imperfect sky-subtraction 
contribute less to the final stacked spectrum. 



2.2 Lick/IDS System 

A lon g term progr a mme un dertaken by Burstein et al.l 
ill984h; I Faber et al. | l|l985l) ; IWorthev et al. I (|1994| ) and 
ITrager et al. 1 1 19981 ). amongst others, has yielded a library of 
stellar spectra, obtained with the Image Dissector Scanner 
(IDS, iRobinson fe Wamplerl I1972T I on the Shane 3m tele- 
scope at the Lick Observatory. These authors have also es- 



tablished a spectral index system used to investigate element 
abundances from low-resolution integrated spectra of extra- 
galactic stellar populations. 

Twenty five absorption features in the wavelength range 
4000 < A < 6000 A at ~ 9 A resolution, have been identi- 
fied that are sensitive to effective temperature (T e g), surface 
gravity (g) and metallicity (Z) of a star. The wavelength def- 
initions of these features were chosen such that the broad ab- 
sorption features found in elliptical galaxies could be studied 
and to minimise the effect of galaxy velocity dispersion on 
the measured line-strengths. 

All features measured on the Lick/IDS system are de- 
fined by a central "feature" bandpass and two adjacent 
"pseudocontinuum" bandpasses, see Appendix[A] From this 
continuum and feature it is possible to measure a "pseudo" 
equivalent width (EW). It is not a true EW as the wave- 
length definitions for the bandpasses are fixed and, depend- 
ing on the element abundances, instrumental resolution or 
velocity dispersion of the stellar population, the wings of an 
absorption feature may extend beyond the feature bandpass. 

As a consequence of adopting fixed wavelengths to de- 
termine the continuum and feature fluxes, the measured line- 
strength depends on the flux in the continua bandpasses 
as well as the feature bandpass. The continua also contain 
absorption features which, when coupled to the effects of 
galaxy velocity dispersions and instrumental resolution, will 
help suppress its average flux. 

Other factors that effect the measured line-strength 
of absorption features are gas, dust and telluric features. 
Although dust has little effect due to the generally nar- 
row nature of the features, gas can be a problem. Neb- 
ular emission will contaminate the measured line-strength 
of the Balmer line s by filling in the "natural" absorption 
l|Moore et alj|2002l . and references therein) and, when using 
the H/3 absorption index, will lead to incorrectly determin- 
ing older ages for stellar populations. While the H/3 index 
line-strength is directly affected by Balmer emission from 
ionized gas, some indices are indirectly affected by emission, 
e.g. the red sideband of Mg b can be contaminated by [Nl] 
iGoudfrooii fc Emsellemlll996r ). 

Telluric emission and absorption can also have a signifi- 
cant impact on the accuracy of the measured line-strengths, 
but stacking spectra in order to improve the S /N of the Lick 
indices also allows us to minimise the impact from these tel- 
luric features. By co-adding spectra at different redshifts, 
telluric features "appear" at different restframe wavelengths 
in the deredshifted spectra. Using the inverse variance of the 
flux when weighting each pixel allows us to minimise the im- 
pact from these features on the final stacked spectrum. 

It should be noted that the SDSS convention is to 
present all wavelengths as vacuum wavelengths. There- 
fore, the wavelength definitions of all Lick Ind ices were 
converted to vacuum wavelengths l|Mortonl [l99lh prior to 
performing any measurements of index line-strengths on 
spectra. We have used the publicly-available software IN- 
DEXF0 l|Cardielll2007l ) to measure our absorption feature 
line-strengths from our co-added stacked LRG spectra. 



http: / /www. ucm.es/info/Astrof/software/indexf/indexf.html 
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Figure 4. Distribution of our quiescent LRG sample in the velocity dispersion-luminosity plane in four different rcdshift slice of Az = 0.1 
over the redshift range 0.0 < z < 0.4. Velocity dispersions have been aperture corrected and absolute magnitudes have been k+e corrected 
to z = 0.2. Contours show the distribution in number density of the whole sample, and are in the intervals 1,10,50,100. The four boxes 
are the same in each panel and show the boundaries used to create the four mass (velocity dispersion) samples used in Section [4] when 
exploring the rcdshift evolution of the SSP parameters for the quiescent sample. 



2.3 Lick Stars in the SDSS 

To properly calibrate to the Lick/IDS system it is necessary 
to have SDSS observations of Lick stars, or galaxies, to allow 
the zero-point offsets to be determined, see Appendix IA1I 
The library of Lick stars that could be used in calibrating 
the SDSS to the Lick system have not been intentionally 
targeted by the SDSS and therefore requires serendipitous 
observations of Lick stars to allow the SDSS to be placed on 
the Lick system. 

The diameter of the SDSS fibres is 3.0" with the median 
seeing at Apache Point of ~ 1.5". We have therefore chosen 
a cut-off of < 1.5" when matching SDSS objects and Lick 
star coordinates. In searching for Lick stars in the SDSS we 
have included both special and survey plates in the search 
and have identified 13 stars from the Lick library that match 
SDSS objects to within 1.5"; 11 stars in M67 (special plate # 
0321) and 2 stars in NGC 7789 (special plate # 2377) and are 
listed in Table [1] The cluster proper motion dispe rsion for 
M67 is ~ 0.8 mas yr" 1 dMcNamara fe Sanderl 197811 and fo r 
NGC 7789 is ~ 0.4 mas vr" 1 l|McNamara fe Solomonlll98ll ). 
with only one star in Table [1] having a measured prope r 
motion of ~ 9.0 mas yr" 1 for M67 F 164 l|Hog et al.lll998h . 
Given the small size of the proper motions we have ignored 
this effect when matching the coordinates of SDSS objects 
to the Lick stars in these two clusters. 

The details of how our SDSS LRG spectra were cal- 
ibrated onto the Lick/IDS system using the Lick stars in 
Table 1 is given in Appendix A. 



3 METHOD 

3.1 Simple Stellar Population Models 

Simple stellar population (SSP) models have become an in- 
valuable tool in allowing the physical properties of stellar 
populations to be probed via measurements of absorption 
features in their integrated spectra. But absorption features 
generally suffer from the same age-metallicity degeneracy 
that effects the interpretation of star formation histories and 
chemical evolution of stellar population using their broad- 



Table 1. Lick stars observed by SDSS 



Star 


Plate 


MJD 


FiberlD 


Spectral Type 


M67 IV- 77 


0321 


51612 


385 


K0 IV 


M67 IV-68 


0321 


51612 


386 


G8 V 


M67 1-17 


0321 


51612 


388 


F0 V 


M67 F 115 


0321 


51612 


463 


F6 


M67 F 105 


0321 


51612 


466 


K2 III 


M67 F 231 


0321 


51612 


479 


K0 III 


M67 11-22 


0321 


51612 


480 


K0 IV 


M67 F 175 


0321 


51612 


490 




M67 F 164 


0321 


51612 


491 


Kl III 


M67 IV-20 


0321 


51612 


499 


K0 III/IV 


M67 F 193 


0321 


51612 


519 


K0 IV 


NGC 7789 676 


2377 


53756 


160 


G8 III 


NGC 7789 897 


2377 


53756 


492 


G9 III 



Note: Spectral classifications obtained from SIMBAD 



band optical colours. Balmer lines become weaker while 
metallic line become stronger as the age and metallicity of 
a stellar population increases. 

While different absorption features behave differently 
to age and metallicity, with some more sensitive to age and 
others to metall icity, a fu r ther complication is that due to 
a-enhancement. I Wort he vl l|l998l ) has shown that regardless 
of isochrones, stellar library, fitting functions or authors, 
the SSP models at the time were unable to explain the 
large spread in Mg at fixed Fe of ellipticals. This overabun- 
dance of Mg is considered to be a consequence of the differ- 
ent timescales involved in the evolution of stars of different 
masses. The abundance of the a-elements (O, Ne, Mg, Si 
amongst others), which are primarily a product of nucle- 
osynthesis in Type II supernovae, is enhanced over the the 
abundance of the Fe-peak elements (Fe, Cr amongst others) 
that result from Type la supernovae. The ages and metal- 
licities derived for SSPs differ depending on whether an a- 
element or Fe-peak element is used to represent metallicity 
([Z/H]). 

Since elliptical galaxies show an enhancement of a- 
element over Fe-element abundance, it has become necessary 
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to employ SSP models that acc ount for these effec ts. In this 
work we use the SSP models of lKorn et all (|2005l , hereafter 
KMT05) which u se the Evo l ution ary Population Synthesis 
(EPS) sch eme of iMarastonl (1 19981 . hereafter M98) and the 
method of iThomas et alj i|2003f ~hereafter TMB03) to pro- 
duce the line-strengths of the Lick indices for stellar pop- 
ulations with variable element abundance. In KMT05, the 
authors use model stellar atmospheres to produce synthetic 
stellar spectra from which they compute line-strengths and 
from this, index response functions as a function of metal- 
licity. They then use the scheme of TMB03 to produce Lick 
index line-strengths that account for variable [a/Fe] ratios. 
These models have been calibrated on globular clusters that 
are known to have enhanced [a/Fe] ratios above the solar 
value. 



3.2 Markov chain Monte Carlo 

Different authors have employed different strategies when 
converting between the observed line -strengths in galax; 



alaxy 
2003; 



spectra and the mod el parameters l|Mehlert et alj _ 
Ide la Rosa et al.ll2007D . However, most previous techniques 
have involved interpolating a grid of model index values at 
the location of the observed index values for the object in 
question. Since the KMT05 models have three SSP param- 
eters, this requires that one of the parameters needs to be 
fixed and there is normally an iterative procedure employed 
to explore the parameter space. At each iteration, a new set 
of model grids are generated and the observed index values 
are used to re-generate the appropriate SSP parameters. In 
order to determine the correct SSP parameters for the ob- 
served index line-strengths, multiple model index grids may 
be used. Moreover, each author may have a different con- 
vergence criteria for these iterations such that the same SSP 
parameters need to be reproduced with different index-index 
grids and are consistent to a pre-determined level. 

In this work, we explore an alternative methodology, 
which is based on a Markov chain Monte Carlo (MCMC) 
technique to perform the iterative minimisation in order to 
obtain the SSP parameters that correspond to the measured 
index line-strengths of the object under study. In the con- 
text of Bayesian inference, the MCMC process produces a 
random sequence of dependent variables that are drawn di- 
rectly from the posterior distribution which is achieved using 
B ayes' Theorem: 



P(%) 



J p(y\0) P (0)d6 



(i) 



where p(0\y) is the posterior probability density and assigns 
a probability to the model, 8, given the data, y, and any 
prior knowledge concerning the model. The probability den- 
sity p(y\8) is associated with obtaining the data given the 
vector of model parameters and is also called the likelihood 
(£). The marginal probability density, p{8), describes the 
probability associated with the parameter vector and en- 
compasses any prior knowledge of the model parameters. 
The denominator in the above expression can be considered 
a normalisation factor and is ignored in our implementation. 

We implement the Metropolis -Hastings algorithm 
^Metropolis et al J 1 19531 : lHastingsl fl97oT) to ensure that the 
stationary distribution of the Markov chain samples from 



the posterior distribution. This form of MCMC uses a can- 
didate generating, or proposal, distribution to propose new 
locations in parameter space which allows the likelihood sur- 
face to be explored. During this exploration the proposed 
new locations are accepted or rejected based on logical cri- 
teria embodied in the Metropolis-Hastings algorithm. 

While the form of the candidate-generating distribu- 
tion docs not impact on the ability of the Metropolis- 
Hastings algorithm to reach a stationary distribution, it does 
imp act on the efficiency o f the convergence to stationar- 
ity (|Doran fc Mulled [20041. The most efficient candidate- 
generating distribution to adopt would be that of the pos- 
terior distribution, but this is not known a priori. In our 
implementation we adopt a multivariate normal distribu- 
tion as our candidate-generating distribution and estimate 
the covariance matrix from the data itself. 

Convergence to the stationary distribution is iden- 

tified through u se of the Gelman-Rubin 7?.-statistic 

|Gelman fc Rubin! Il992p with convergence identified when 
1Z ^ 1.1 (|Verde et alJl2003T ) for each parameter of the model. 
Identifying convergence allows the "burn-in" phase to be 
excluded from the chain and accurate estimation of the 
confidence-intervals on the model parameters. 



3.3 Likelihood Estimation 

In order to use the SSP models of KMT05 with our MCMC 
approach we adopt a trilinear interpolation scheme, since 
the models have three parameters (age, metallicity and a- 
iron ratio), such that the SSP models can be investigated 
at an arbitrary location in parameter space. While there is 
a unique mapping from model parameters to index values, 
there is no unique mapping from index values to model pa- 
rameters and therefore the models cannot be inverted. We 
must therefore use multiple Lick indices simultaneously in 
order to break this degeneracy. 

The Metropolis-Hastings algorithm employs a likeli- 
hood ratio test in order to determine acceptance or not of 
a candidate position in parameter space. The likelihood 
(C = p(y\8)) is obtained by assuming that the data, y, 
is a set of N independent normally dis tributed random 
variables, yt,i = 1,...,N l|Cowanl 1 19981 ). Each yi has a 
different and unknown mean, fii, but known variance, 
of. The joint p.d.f. is the product of these N normal 
distributions: 



iV 



exp 



-(y» - mf 
2a? 



(2) 



where the true values, fu, depend on the model parameters, 
i.e. 8 = . . . , 9 m ). Taking the logarithm of the joint p.d.f. 
yields 



m£(0)=V — 



i^(vi-fMi(0)Y 



(3) 



We can maximise the log-likelihood by minimizing the \ 2 
function, 
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N 



( yi - M *(0)) 2 



(4) 



In the context of estimating SSP parameters from line- 
strength data, yi ±<7j = Ii ± d L, where L is the line-strength 
of a particular absorption feature, dl« is its one-sigma error 
and i represents the Lick index used, e.g. H/3 , Mgb etc. 
The "true" absorption indices, are determined from in- 
terpolating the SSP models of KMT05 for a particular set 
of model parameters, = {t, [Z/H] , [a/Fc] }. 

At each location in parameter space, we use the SSP 
models to determine the corresponding values for the cho- 
sen Lick indices using the the trilinear interpolation method 
to go from model parameters to Lick index line-strengths. 
The MCMC approach is then used to iteratively maximise 
the likelihood, such that at each stage in the evolution of 
the chain, the likelihood, £ n +i, calculated at the new can- 
didate position, n +i, is compared with the likelihood, C n , 
determined at the present location in parameter space, n . 
The Metropolis-Hastings algorithm is then used to make the 
decision on accepting or rejecting the candidate position. If 
the next step is rejected then the current step is re-saved as 
part of the chain. 

The chain now performs a random walk in parame- 
ter space and generates the sequence of parameter samples 
{0i, 02, . . . , On, 0„+i}. Once the chain has converged to the 
stationary distribution, i.e. burned in, the random walk will 
be confined to the vicinity of the global mode in the posterior 
distribution. In Appendix [Bl we test this method of SSP pa- 
rameter esti mation by comparing with previously published 
results from lThomas et alj (|2005l . hereafter TMBO). 



4 RESULTS 

4.1 Line Strengths 

The relationship between the line-strength and redshift for 
the four different mass (velocity dispersion) samples, defined 
in Figure [4] are presented in Figure [5] for a selection of the 
key Lick indices used in age-dating of galaxy spectra (again 
see Appendix for details of how the SDSS spectra were cali- 
brated onto the Lick system) . These relationships are consis- 
tent with that of a passively evolving stellar population, i.e. 
the strength of H/3 and HyF indices increase with redshift 
while the Mg b and (Fe) lines decrease in strength, indicated 
by the dotted line in each panel. 

The scatter in these relationships is consistent with 
the error on a typical line-strength and there is a clear 
trend with velocity dispersion, i.e. Mgb and (Fe) show 
a positive correlation with velocity dispersion while H7J 7, 
shows a negative correlation. These trends are consistent 
with t he finding from o t her a u thors, e.g. | Bender et al.l 
J1993T): iKuntschner et all (|200ll ); iBernardi et alj (|2003bl ); 
lLoubser et al.l ( 20091 ). While there is no obvious segrega- 
tion of the H/3-redshift relationships there is a hint that 
the 260 < a v sC 290 km s _1 (circles/ green) sample yields 
line-strengths that are systematically larger than the 320 < 
a v ^ 350 kms _1 (diamonds/pink) sample, consistent with 
the negative correlation observed by other studies. 



4.2 SSP Parameter Estimates - Uniform Priors 

Using the H7F, Mgb and (Fe) data in Figure [5] and 
the parameter estimation method outlined in Section [3] 
we determine the SSP parameters for each mass sample 
and the results are shown in Figure These parameter 
estimates have been obtained using uniform priors that 
cover the full parameter space of the KMT05 models, i.e. 
age (Gyr) ~W(0-31, 18.2), [Z/H] ~ U(-1.0, 0.97), [a/Fe] ~ 
U{— 0.25, 0.73), when computing the posterior probability 
using Equation [TJ and then estimating the 16, 50 and 84 
percentiles from the posterior distribution for each param- 
eter after marginalising over the remaining parameters. In 
the rest of this paper, the median is used to define the "lo- 
cation" and the 16 and 84 percentiles the la error on SSP 
parameters. 

The evolution in [a/Fe] with redshift is somewhat con- 
fusing as the low (blue triangles) and high (pink diamonds) 
mass (velocity dispersion) samples exhibiting significant 
scatter and may not be consistent with a constant value 
over the entire redshift range (there is some indication for a 
change in behaviour at z > 0.35). The intermediate mass 
samples (green circles and red squares) are more stable 
and exhibit no significant gradient, as expected for galax- 
ies that are not chemically evolving. There does however 
appear to be a positive correlation between velocity disper- 
sion and [a/F e ] whic h is consistent with other studies, e.g. 
iThomas et all (|2005h . 

The evolution in [Z/H] with redshift is also consistent 
with a chemically unevolving population of objects, al- 
though the lowest mass sample does exhibit a significant 
gradient. The metallicity also exhibits a significant positive 
correlation with velocity dispersion, where the more massive 
objects have higher metallicity, again consistent with other 
studies. Finally, the evolution in age with redshift exhibits 
the expected behavior; galaxies become younger with red- 
shift, with a trend that suggests the most massive objects 
are, in fact, the youngest. 

While the H/3 and H7.F 1 trends are consistent with trac- 
ing a passively evolving population, it is not possible to ob- 
tain consistent estimates of the SSP parameters using these 
indices. The dotted line in Figure [S] is the relationship ex- 
pected for a population of objects in a ACDM concordance 
cosmology, with [Z/H] = 0.37 and [a/Fe] = 0.27 and a for- 
mation age of 4.5 Gyr. To reproduce the H/3 line-strengths, a 
population with the same chemical composition would need 
a formation age of approximately 2.5 Gyr and hence yield 
inconsistent line-strengths for H7. 

In Figure [7] we show the correlations between the pa- 
rameter estimates of Figure [S] and in Table f2] we show the 
Spearman rank correlation coefficient ( — 1 ^ p ^ +1) and 
in brackets the significance of the correlation from zero 
(0 ^ s p ^ 1; with low values having greater significance). 
From this we can see that, except for the lowest mass sam- 
ple, correlations between the parameters have little statisti- 
cal significance, i.e. values for s p > 0.1. While interpretation 
of the lowest mass sample is difficult because of the gradients 
in [Z/H] and [a/Fe] with redshift, possibly due to including 
a population of young objects at low redshift (c.f. the H/3 
redshift relationship in Figure [5}, there does appear to be a 
significant correlation between the chemical composition of 
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Figure 5. Line-strengths of H/3, Mgb , (Fe) and H7F as a function of redshift for the four samples of Figure 3; 230 < a v ^ 260 km s 1 
(blue triangles), 260 < a v sC 290 km s" 1 ( green circles), 290 < cr v ^ 320 km s 1 (red squares), 320 < a v ^ 350 km s 1 (pink diamonds). 
The dotted line shows the expected variation in the line-strengths for an object with [Z/H] = 0.37 and [a/Fe] = 0.27 and a formation 
age of 4.5 Gyr for a ACDM cosmology. Segregation of the line-strengths for the different samples is obvious except in the case of H/3. 
The more massive the sample (i.e. the higher velocity dispersion), the stronger the Mgb and (Fe) line-strengths and the weaker the H7F 
line-strength at a given redshift. Errors are only provided on the the 260 < cr v ^ 290 kms - 1 sample (green circles) to avoid overcrowding 
on the plot. The errors on the other relations are similar as they have comparable signal-to-noise. 



Table 2. SSP parameter correlations 



Sample 




Total Number 


t-[Z/H] 


H 


01 /Fe] 


[a/Fe]- 


-[Z/H] 


230 < <j v 


C 260 kms -1 


4472 


0.54 (0.05) 


-0.71 


(0.01) 


-0.64 


(0.02) 


260 < a v 


^ 290 kms" 1 


6195 


-0.26 (0.33) 


-0.15 


(0.58) 


0.20 


(0.45) 


290 < a v 


^ 320 kms- 1 


4835 


-0.17 (0.57) 


0.05 


(0.85) 


0.24 


(0.41) 


320 < a v 


^ 350 krns" 1 


2350 


-0.55 (0.16) 


-0.02 


(0.96) 


0.36 


(0.39) 



these objects and their mass, the more massive the objects 
the higher their [Z/H] and [a/Fe] . 

In Figures 16171 we have employed H7F, Mgb and 
(Fe) when estimating the SSP parameters. In Figure [S] we 
show the results of using H/3 instead of H7F when re- 
constructing the age-redshift relationship. The [Z/H] and 



[a/Fe] redshift relationships (not shown) are similar to those 
in Figure [6] but shifted to lower values by approximately 
0.05 dex in [Z/H] and 0.01 dex in [a/Fe] . The scatter in R/3 
translates into significant scatter in the age estimates but 
the same trend seen in Figure [6] still exists, i.e. the age of 
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Figure 6. Reconstruction the evolution of the mean luminosity-weighted age, [Z/H]and [a/Fe] using H7F, Mgb , Fe5270 and Fe5335. 
Symbols and colours are the same as shown in Figure [5] Open symbols represent data with S/N< 100 A -1 in H7F. Horizontal lines 
represent the mean [Z/H] and [a/Fe] for data with S/N> 100A _1 in the 230 < a v ^ 260 kms" 1 (blue short dash), 260 < a v < 290 kms" 1 
(green dot-dash), 290 < cr v ^ 320kms _1 (red dotted-dash) and 320 < a v ^ 350kms — 1 (pink long dash) samples. For the reconstruction 
of the age-redshift relationship the dotted line shows the age of the universe t\j(z) for a ACDM cosmology while the dashed line indicates 
tjj(z) — 4.5 Gyr and is for reference only. To reduce clutter in the diagram we again only show the la errors on the 260 < <r„ ^ 290kms — 1 
(circles/green) sample only, but errors for each sample are consistent with the scatter. 
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Figure 7. Correlations between the SSP parameters of Figure [6] Parameter estimates are formed from the one-dimensional posterior 
distribution marginalised over the remaining parameters. 




0.4 



Figure 8. Reconstruction the evolution of the mean luminosity- weighted age using H/3, Mgb , Fe5270 and Fe5335. The only difference 
between this figure and Figure [6] is the choice of Balmer line used in the estimation of the SSP parameters. While the trend of the age- 
redshift relationship is the same as that derived using H7F, there is more variability in the age estimates resulting from the variability 
in H/3 as seen in Figure[5] Also the ages derived using H/3 are systematically older than those derived using H7F. Symbols and lines are 
the same as in Figure \E\ and we have excluded the [Z/H] and [a/Fe] traces here as they are similar to the corresponding traces in Figure \E\ 



the objects decreases with redshift, although the ages are 
older by approximately 2.5 Gyr. 

4.3 SSP Parameter Estimates - Gaussian Priors 

By binning our samples into redshift intervals, we have ob- 
tained multiple independent estimates of the chemistry of 
our LRG population. For each mass sample we can use 
these independent estimates of [Z/H] and [a/Fe] to produce 
a prior probability distribution, which we can then employ 
to better constrain the age estimates when estimating the 
posterior probability using Equation [T] When determining 
the priors on [Z/H] and [a/Fe] , we only use individual es- 
timates that have been produced using line-strength data 
with S/N> 100 A" 1 in H7F 1 and we generate Gaussian 
priors from the mean and standard error (rms/^/n) of the 
[Z/H] and [a/Fe] distribu tions for a given mas s sample. We 
use importance sampling (|Lewis fc Bridlef 2002) and the new 



Gaussian priors to re-weight the original MCMC chain data 
obtained with uniform priors on the parameters. Although 
[a/Fe] does not correlate with age, [Z/H] does (see Figure 3 
of TMBO) and placing a prior constraint on [Z/H] results in 
a tightening of the multivariate posterior distribution yield- 
ing a more localised marginalised distribution for the age of 
the object. 

In Figure [9] we show what effect the Gaussian priors 
on [Z/H] and [a/Fe] have on the reconstruction of the age- 
redshift relationship for each mass sample. The impact of the 
prior on [a/Fe] is negligible, but the decrease in the errors on 
the individual age estimates, and the reduction in the scatter 
for a given mass sample resulting from the prior on [Z/H] , 
and the age-metallicity degeneracy, is significant. Over the 
redshift range studied herein, we see a decrease in the age of 
LRGs of ~ 5Gyrs, which is fully consistent with expectations 
from a ACDM universe. 
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Figure 9. Reconstruction the evolution of the mean luminosity- weighted age using H7F, Mgb , Fe5270 and Fe5335 and employing 
Gaussian priors on [Z/H] and [a/Fe] . The size of the error bars on the 260 < a v ^ 290 kms -1 (circles/green) sample are about the same 
size as the data points. 



5 DISCUSSION 

We present in Figures 4, 5, 7, 8 and 9, the results of our 
analysis of LRG ages, metallicities and a-enhancements as 
a function of redshift. We also provide in Appendices C and 
D the data presented in these figures. Overall, these trends 
are consistent with expectations and can provide important 
constraints on cosmology (via the age-redshift relation) and 
galaxy evolution studies. However, we raise here a number 
of important caveats that should be considered when using 
these data. 

First, we have tried to select a consistent population 
of quiescent galaxies with redshift using the SDSS Lumi- 
nous Red Galaxy (LRG) selection, which is designed to se- 
lect "red and dead" massive galaxies out to z ~ 0.4 (with 
a constant number density). The size of this LRG sample 
also allows us to preferentially select quiescent galaxies with 
no obvious emission lines, as well as define four subsamples 
with the same range of velocity dispersions and luminosities. 
However, we have employed a simple evolutionary correction 
that maybe too naive at higher redshift and also performed 
aperture corrections to velocity dispersion measurements for 
which it has been necessary to assume that corrections de- 
rived at low redshift objects can be applied to higher redshift 
objects. This assumes no evolution in velocity dispersion, or 
stellar population distributions, and hence no evolution in 
the mass distribution, over the range 0.0 < z < 0.4. Both 
of these assumptions may impact on the samples used in 
generating the stacked spectrum and could potentially re- 
sult in different populations at high and low redshift being 
probed. This could explain the gradients in the [Z /H] and 
[a/Fe] -redshift relationships even if each population probed 
is passively evolving. 

While some studies into the luminosity function for 
the red sequence have shown that the most massive sys- 
tems show little or no evolution over the redshift rang e 
probed in this work, (See lBell et alJl2004lFaber et al.ll2007f ). 
this may not be true for lower mass systems in our sam- 



ple. The anticorrelation between [Z/H] and [a/Fe] for the 
230 < cr v ^ 260 sample suggests that the low mass sys- 
tems become less alpha-enhanced and more metal-enriched 
at late times. This could arise if our galaxy sample includes 
a fraction of "blue cloud" galaxies whose star formation is 
switching off and hence moving onto the red sequence since 
z=0.4. These low mass systems have been forming stars over 
longer periods of time and have experienced more chemical 
enrichment than galaxies with similar velocity dispersions 
whose star formation switched off at earlier times, giving 
them lower [a/Fe] and higher [Z/H] at low redshift. This low 
mass sample, while possibly offering insights into the mech- 
anisms behind galaxy evolution, is less useful in its intended 
role as a cosmological probe due to the possibility of not 
accurately tracing a passively evolving sample of galaxies. 

There are further complications from adopting the 
Lick/IDS system when estimating the SSP parameters, even 
though this is a popular and well-known method. For exam- 
ple, matching the SDSS instrumental resolution to that of 
the IDS instrument, as well as performing each of the aper- 
ture, velocity dispersion and zero-point corrections, poten- 
tially introduces a source of error into our line-strengths and 
thus SSP parameter estimations. Also, the KMT05 models 
are not perfectly matched to the Lick/IDS system and there- 
fore, the accuracy of our SSP parameter estimates could de- 
pend on the difference between the calibration of the data to 
the models. This difference varies for each of the Lick indices, 
resulting in the absolute value of any SSP parameter being 
dependent on the particular set of Lick indices chosen. This 
is illustrated in Figures [6] & [8] by the discrepancy between 
the absolute ages for our quiescent LRGs when derived us- 
ing H7 or H/3 . Therefore, we would caution the reader from 
using our absolute ages in this paper, but recommend they 
focus on the relative ages with redshift. 

We also note that our observed line-strengths are 
strongly dependent on the accuracy of the velocity disper- 
sion correction (see Appendix A). Since the velocity disper- 
sion correction depends on the stellar sample used to gener- 
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ate the correction factors, matching the stellar sample to the 
dominant stellar population of the LRG may become impor- 
tant in the future. Also, at higher redshift, the contribution 
of younger stellar populations to the galaxy spectrum may 
require a different stellar sample to be used in making the 
velocity dispersion correction. Therefore, questions remain 
over the appropriateness of employing a single stellar sample 
of Kill stars to determine the velocity dispersion correction 
over the redshift range used in this work. 

However, we finish by stressing that for our samples 
we see little evolution in the shape of [Z/H] and [a/Fe] with 
redshift, especially for the intermediate mass (velocity dis- 
persion) samples which are consistent with a constant over 
the redshifts probed. This is re-assuring as it does indicate 
that our attempts to select a passively evolved, quiescent 
subsample of galaxies that could be used to probe the age- 
redshift relation can be found and the errors well-controlled. 
We provide these relationships here for others to use in their 
analyses, and plan to explore the cosmological constraints of 
our age-redshift relationships in a subsequent paper. 



6 SUMMARY 

We present here a detailed analysis of a sample of passive 
SDSS LRGs which are selected to be the same population 
of galaxies over the redshift range 0.0 < z < 0.4. A total 
of 17,853 LRGs are co-added in four bins of velocity dis- 
persion and luminosity to provide high signal-to-noise spec- 
tra for detailed spectra absorption line measurements. In 
Section 3 and Appendix A, we outline the careful calibra- 
tion of these absorption line measurements onto the well- 
established Lick/IDS system which has been extensively 
used in the literature to constrain Simple Stellar Population 
(SSP) models and thus derive galaxy parameters like age, 
metallicity and Q-enhancements. In Figures 4, 5, 7 and 8, 
we present our results which show clear trends for these pa- 
rameters with redshift (we also provide the data from these 
figures in Appendices C & D). In particular, for our two 
intermediate mass (velocity dispersion) samples, we see a 
constant [Z/H] and [a/Fe] with redshift which confirms our 
assumptions that these LRGs are a passively evolving sub- 
sample of galaxies (i.e., they are chemically unevolved over 
this redshift range) providing confidence on our age determi- 
nations for these galaxies. We see a clear trend of decreasing 
age of our LRGs as a function of redshift (~ 5Gyrs over the 
redshift range probed here), which is fully consistent with 
expectations from a ACDM universe. It also appears that the 
most massive sample of LRGs is also the youngest (Figures 
5 and 8). We provide these relationships now to help others 
studying cosmology and galaxy evolution, as well as provide 
in our appendices, the methodology required for others to 
calibrate SDSS spectra onto the Lick/IDS system. 
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APPENDIX A: 
SYSTEM 



CALIBRATING TO LICK/IDS 



We provide here the details of our calibration of SDSS spec- 
tra onto the Lick/IDS system. We use the most recen t def- 
initions of the Lick indices from Tra ger et al.l (|l998T ) and 
IWorthev fc Ottavianil (|l997T ). while the index line-strengths 
are calculated using the standard equations 



^ (A) 
and 

Jm (mag) 



/" A c 2 




/a c1 


C(X) 



dA. 



-2.51o gl 



J^[F(X)/C(\)]dX- 
Ac 2 — A Cl 



(Al) 



(A2) 



where the narrow feature atomic indices, /a, have their line- 
strengths expressed in Angstroms (A) while the broad fea- 
ture molecular indices, Im, are expressed in magnitudes 
(mag). Also, in Equations IA1I and IA2I \c 1 and Ac 2 are the 
wavelength definitions of the feature bandpass and F(A) and 
C(A) are the flux in the index bandpass and pseudocon- 
tinuum respectively. When measuring the line-strengths of 
absorption fea tures we use the software package INDEXF 
l|Cardielll2007l ). 

The flux in the local pseudocontinuum, C(A), is deter- 
mined using 



1 Ar — Ai Ar — Au ' 



where 



X Xh -' Xb - 



and 



Ft s f 2 F W d A. 

and 
A 



b = ( A b 2 + A b 1 )/ 2 ' A r = (Ar 2 +A ri )/2. 



(A3) 



(A4) 



(A5) 



(A6) 



A bi , A v 



, \r ± and Ar 2 are the wavelength definitions of the 
blue and red bandpasses. The resulting local pseudocontin- 
uum, C(A), generated by this approach is a straight line 
that connects the average flux, defined at the central wave- 
length, of the blue (Fj^Ak) and red (Fr,Ar) bandpasses of 
the index. 
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Al Transforming to Lick/IDS 

Since the Lick indices have been defined with spectra taken 
on the IDS, the absorption index system developed by the 
Lick group has the characteristics of the IDS instrument in- 
herent to it. Also, the spectra obtained by the Lick group 
were not flux calibrated and this results in the shape of 
the continuum of a Lick/IDS spectra differing to the true 
shape. Line-strengths measured on spectra not observed 
with the IDS instrument need to be transformed to the 
Lick/IDS system such that they can be compared with 
Lick/IDS calibrated SSP models. For galaxy spectra this re- 
quires that spectra are smoothed to the IDS resolution and 
line-strengths are corrected for aperture and velocity disper- 
sion effects and zero-point corrections applied. These steps 
are described in more detail in the following subsections. 

A2 Matching Instrumental Resolutions 

The measured line-strength of a given Lick Index depends 
on, amongst other things, the instrumental resolution. As a 
result of the fixed wavelength definitions of the index band- 
passes and pseudocontinua, the line-strength generally de- 
creases as the instrumental resolution is degraded. This in- 
strumental broadening causes the wings of the spectral fea- 
ture to extend outside their bandpass window reducing the 
strength of the absorption feature within the window. In 
addition to this, adjacent spectral features may depress the 
flux in the continuum by being broadened into the contin- 
uum windows reducing the line-strength further. 

In order to transform SDSS spectra to the Lick/IDS 
system, it is necessary to degrade the higher resolution SDSS 
spectra to match the lower resolution of the IDS instrument. 
But, because the Lick/IDS calibrated SSP models produce 
index strengths for a stellar population in the rest frame 
at the resolution of the Lick/IDS system, it is necessary to 
ensure that the absorption feature in the observed spectrum 
has the same instrumental resolution as it would have at 
rest in the Lick/IDS system. 

The SDSS spectrographs have 640 fibers with the spec- 
tral resolution of each fiber having a slightly different depen- 
dence with wavelength. The response of the IDS instrument 
can be treated as a Gaussian of wavelength dependent width 
jWorthev fc Ottavianill 19971 ) and to match instrumental res- 
olutions, the SDSS spectra were convolved with a Gaussian 
whose wavelength dependent width was found using 

OLickW 1 + z ) ^Lick) = ""SDSsOobs) + ^match (\>bs) , ( A 7) 

where cr ma tch(^obs) ls the width of the Gaussian required 
to match the Lick/IDS and SDSS spectral resolutions, 
°"Lick((l + z ) ^Lick) is the wavelength dependent spectral res- 
olution of the Lick/IDS instrument shifted to the observed 
frame and crgrjss (-^observed) is the SDSS resolution for a par- 
ticular fiber. 

The SDSS spectral resolution is determined for each 
object separately using information in the header (HDU #6) 
of the FITS file for a given object. This header contains the 
RMS resolution (in pixels) as a function of pixel number 
and to convert from RMS resolution in pixel units to RMS 
resolution in wavelength units, the RMS resolution in pixels 
is multiplied by the local pixel size in wavelength units 

°\ = CTpbc x ln(10) x A pix x 10" 4 , (A8) 



where A p ; x is the wavelength at a given pixel and the fac- 
tor 10 -4 is a result of the SDSS spectra being binned with 
constant logarithmic dispersion of this size. 

The convolution of the SDSS spectra with a Gaussian 
was performed in pixel space and is given by 

TO — 1 

(S ® G{j)) 3 = S j+i - m/2 Gi(j) , (A9) 

i=0 

where S is the science spectrum, G(j) is the variable width 
Gaussian and m is the width of the Gaussian window func- 
tion and is set at 4.5 times the FWHM of the Gaussian. 
The width of the Gaussian is wavelength dependent, which 
in pixel space means that the width is a function of pixel 
location j. At each pixel j, a normalized Gaussian is gener- 
ated whose width in pixels is determined from Equations lA7l 
and lA8l where £>"Lick(-M is inferred from linear in t erpo lation 
of the data in Table 8 of IWorthev fc Ottavianil l|l997h . af- 
ter converting the air wavelengths in the table to vacuum 
wavelengths and redshifting to the frame of the observed 
spectrum, A obs = (1 + z) A Lick . 

As a result of this convolution procedure, the resolution 
of the SDSS spectra at the observed wavelength is matched 
to the resolution at the corresponding rest-frame wavelength 
on the Lick/IDS system, and consequently absorption fea- 
tures measured in the target spectrum match the resolution 
of the same feature on the Lick/IDS system. 

A3 Velocity Dispersion Corrections 

The observed spectrum of a galaxy is equivalent to the con- 
volution of the instrumental response function, the distri- 
bution of the line-of-sight velocities of its stars and their 
integrated spectrum. The velocity dispersion of the galaxy 
and the instrumental response both broaden spectral fea- 
tures, thereby depressing the measured line-strengths from 
their intrinsic value. 

The Lick/IDS system is based upon an empirical stel- 
lar library and is therefore defined at zero intrinsic velocity 
dispersion. In order to compare galaxy index measurements 
to their SSP model predictions, the raw measurements need 
to be corrected to zero velocity dispersion. 

To ge nerate the corre c tion f actors, we adopted the ap- 
proach of IPoggianti et all (|200ll ). which involves matching 
the resolution of a sample of SDSS stellar spectra to the 
Lick/IDS resolution. Then, each stellar spectrum was fur- 
ther convolved with a Gaussian, in logarithmic wavelength 
space, to mimic the velocity dispersion of a galaxy and 
was repeated for velocity dispersions in the range a = 
- 420 km s" 1 in steps of 20 km s" 1 . 

When estimating velocity dispersions, the SDSS 
pipeline uses template spectra which are convolved to a max- 
imum velocity dispersion of 420kms~ 1 , which we then adopt 
as the maximum velocity dispersion to broaden our stellar 
spectra by. Incidentally, the template spectra used by the 
SDSS pipeline are stellar spectra taken from plate 321, the 
same plate used in this work to calibrate the systematic off- 
sets and place the SDSS on the Lick/IDS system. 

The selection of the stellar spectra used to generate 
the correction factors utilises the SEGUE Stellar Parame- 
ter Pipeline (SSPP; lLee et al.ll2007T l database via the CAS. 
As part of the SSPP, the stellar spectra are cross-correlated 
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Table Al. Polynomial fits to the velocity dispersion correction 
factors 



Index 


bo 






&3 






(xl0~ 3 ) 


(xl0~ 6 ) 


(xlO -9 ) 


CN1 


0.000 


0.002 


0.077 


-0.071 


CN2 


-0.000 


0.002 


0.226 


-0.222 


Ca4227 


0.991 


0.420 


0.259 


11.978 


G4300 


0.999 


0.034 


0.390 


0.220 


Fe4383 


0.998 


0.035 


2.035 


-0.803 


Ca4455 


0.986 


0.249 


2.919 


3.582 


Fe4531 


0.998 


0.039 


1.401 


0.526 


C4668 


1.000 


-0.013 


1.148 


-0.037 


Hbeta 


1.003 


0.037 


-0.256 


1.105 


Fe5015 


0.999 


0.075 


2.682 


-1.948 


Mgl 


0.000 


-0.001 


0.076 


-0.091 


Mg2 


-0.000 


0.004 


0.020 


0.029 


Mgb5177 


1.001 


-0.041 


2.338 


-0.555 


Fe5270 


1.001 


0.064 


2.687 


-1.730 


Fe5335 


1.004 


0.033 


4.967 


0.823 


Fe5406 


1.004 


-0.074 


5.100 


1.181 


Fe5709 


0.995 


0.146 


0.905 


3.824 


Fe5782 


0.999 


0.105 


3.533 


4.859 


Na5895 


1.002 


0.051 


0.259 


2.013 


TiOl 


-0.000 


-0.001 


0.023 


-0.021 


Ti02 


-0.000 


-0.000 


0.016 


-0.013 


HdA 


0.999 


0.013 


0.996 


-0.421 


HgA 


1.000 


-0.013 


-0.081 


0.918 


HdF 


0.991 


0.054 


-1.457 


1.398 


HgF 


1.001 


0.019 


0.958 


-1.437 



with ELODIE () high resolution spectra, degraded to match 
the SDSS resolution, to determine radial velocities. Those 
spectra whose best matching template was a K0-K3 III star 
were used to create a sample from which the velocity dis- 
persion correction factors were produced. 

The index strengths for each star at each broadened ve- 
locity dispersion were measured and compared to the index 
strength from the zero-velocity dispersion stellar spectra. 
The correction factors Ci(a) — Index((j = 0)/Index(<r) for 
atomic indices, and Ci(a) — Index(a = 0) — Index(cr) for 
molecular indices, were generated from the broadened spec- 
tra. 

We then fit a polynomial of the form 

ci{o)= J2 h ° l - ( Al °) 

^ i ^ 3 

to the correction factor data for each index where a repre- 
sents the velocity dispersion. This polynomial fit is employed 
when making velocity dispersion corrections to a given line- 
strength. In Figure |A"T1 we show the correction factors and 
the best fitting polynomial for each index. In Table IA1I we 
present the coefficients of the best fitting polynomial to Ci 
for each Lick index. 

A4 Zero-Point Corrections 

The Lick/IDS spectra were not flux calibrated but in- 
stead normalized to a calibration lamp. Differences between 
the continuum of the fluxed SDSS spectra and normalized 
Lick/IDS spectra will manifest themselves as systematic off- 
sets between our measured line strengths and those of com- 
mon objects already on the Lick/IDS system. To determine 



Table A2. Lick index zero-point offsets 
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these offsets, and hence facilitate the use of SSP models cal- 
ibrated to the Lick/IDS system, it is necessary to measure 
the line strengths of Lick stars that are common to SDSS. 

Using spectra of the 11 stars in M67 and 2 stars in NGC 
7789 previously identified as Lick stars, we have matched res- 
olutions to the Lick/IDS and measured line-strengths for all 
25 Lick indices. In Figure IA2I we sh ow the comparison be - 
tween our measurements and those of IWorthev et all (|l994T ) 
for the 13 Lick stars observed by SDSS. In Table EH the 
statistical significance of the mean offset (| offset\/a) for the 
majority of indices can be seen to be less than 2a, where 
a = rms/\/~N. Consequently, we choose only to apply zero- 
point offsets to those indices with offsets that have > 2a 
significance and these are shown in Figure lA"2l as panels with 
bold borders. 



A5 Aperture Corrections 

The measured velocity dispersion and the line-strengths of 
absorption features for a galaxy is affected by a number 
of factors which include the size and shape of the spectro- 
graph aperture used to sample the light from the galaxy, the 
distance to the galaxy, the luminosity distribution and the 
distribution of stellar velocities. 

Elliptical galaxies exhibit radial gradients in their ve- 
locity dispersions as well as in their absorption-line indices 
l|Mehlert et al-l fcoOO), and it is necessary to correct for these 
gradients when comparing our results to other results in the 
literature. Corrections are also necessary when comparing 
data over a large range in redshift as the fixed spectrograph 
aperture will sample different physical scales depending on 
the distance to the galaxy. 

A further complication when comparing results to those 
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Figure Al. Velocity dispersion correction curves for each Lick index, determined from a sample of 31 K III stars. The size of the symbol 
at a given a indicates the magnitude of the standard deviation in correction factors for the stellar sample used. 



in the literature is the size and sha pe of the spectrograph 
aperture used. We use the method of lJorgensen et all (|l995l ) 
to convert rectangular apertures to an equivalent circular 
aperture of radius r ap - They show that a velocity dispersion 
estimated with a circular aperture would yield the same 
mean value, to within 4%, as that through a rectangular 
aperture if 2r ap = 1.025 x 2(xy /tv) 1 ^ 2 , where x and y are the 
dim ensions of the re ctangular aperture. 

|j0rgensen et al.l also show how the velocity dispersion, 
measured with a circular aperture, varies with aperture ra- 
dius and they have determined that within an effective ra- 
dius, the velocity dispersion profile normalised to <r e g, can 
be approximated by a power law of the form 

<A11) 

where a is the gradient in the velocity dispersion, r e is the 
effective radius of the galaxy and <r e g is the velocity disper- 
sion measured through an aperture of r e /8. 

Using Equation I Al II and following the approach of 
iBernardi et al.l |2003af ). we correct the measured velocity 
dispersion, a ap , returned by the SDSS pipeline to a veloc- 



ity dispersion, a corr = cr e 8, that would be measured through 
r e /8. From r e = {b/ a^^rdev , we determined r e , where rdev 
is the de Vaucouleurs radius, b is the semi-minor axis and a 
is the semi-major axis of the galaxy, i.e. r e is the effective cir- 
cular radius. We use the values of rdev, a and b determined 
from the SDSS r-band photometry and obtained using the 
parameters deVrad_r and deVAB_r from the CAS. 

Similar expressions to Equation lAlll are used to correct 
the measured line-strengths to the standard-size aperture 

log I CO rr = log I Tap ~ &I log — ^ , (A12) 

for atomic indices and 

Icorr = Ir ap ~ a I log — ^j- , (A 13) 

for molecular indices, where ai is the radial gradient 
(A log 1/ A log r) for the index and I Tap is the line-strength 
measured from the SDSS spectra. 

In order to correct the velocity dispersions and index 
measu rements for aperture effects we adopt the index gradi- 
ents of lMehlert et all (|2003T ). in which they have determined 
gradients for H^, Mgb, Fe5270, Fe5335 and the velocity dis- 
persion, a, for a sample of 35 early-type galaxies. 
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Figure A2. Comparison between our measurements for the 25 indices of Lick/IDS system and those of IWorthev et a~ | (QUI), for 11 
stars in the open cluster M67 (closed symbols) and 2 stars on the open cluster NGC 7789 (open symbols) observed by SDSS. The mean 
offset for each index is indicated by the solid line. Some panels have less that 13 data points because of missing data from either the 
IWorthev et al.l dataset or our own. Indices whose systematic offset is determined to be > 2 a have had the correction applied and are 
shown as having bold borders. The median error bar is indicated for each index and A Index is in the sense SDSS - Worthcy. 



A6 Comparison with Previous Data 

In order to test our transformation/calibration to the 
Lick/IDS system it is necessary to perform similar compar- 
isons to that already used in determining the zero-point off- 
sets. There exists extensive work in the literature where au- 
thors have obtained line-strengths determined from galaxy 
spectra and have also calibrated their line-strengths to the 
Lick/IDS system. Not every source in the literature explores 
the full set of Lick indices so we restrict our comparison to 
LL3 , Mgb , Fe5270 and Fe5335. All objects that have been 
selected for study in this comparison have been matched to 
within 5" of an SDSS object and have a velocity dispersion 
available from the SDSS pipeline. 

M oore et all (|2002l ) have observed early-type galaxies in 
Coma of which 49 have SDSS spectra. All data, both liter- 
ature and SDSS, have been corrected to a standard circular 
aperture of r e /8, where r e has been determined for each 
object from their SDSS r-band photometric parameters and 
the method outlined in Section lA"5l From Figure [A"3l and Ta- 



Table A3. Testing the calibration to Lick/IDS System using 
Coma galaxies common to SDSS and Moore et al. 



Lick Index 


Objects with S/N > 30 A - 1 


All Objects 


H/3 


+0.05 


+0.06 


Mgb 


+0.03 


-0.06 


Fe5270 


+0.01 


+0.02 


Fe5335 


+0.02 


-0.04 



ble IA31 it would seem that our calibration to the Lick/IDS 

system is within approximately 0.05 A. 

Given that the 2.7" fiber used byH oore et all matches 
well to the 3" fiber used by SDSS, aperture corrections 
between these two datasets will be a minimum. The ve- 
locity dispersion corrections derived by iMoore et alj also 
match to those derived in this work to better than ~ 1% 
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Figure A3. Comparison of line-strengths for H/3 , Mgb , Fe5270 and Fe5335 with iMoore et all (120021 ) for 49 early-type galaxies in 
the Coma cluster. All SDSS spectra have been matched to the Lick/IDS resolution, line-strengths corrected (both SDSS and literature 
measurements) to a standard circular aperture of r e /8 and have been corrected for velocity dispersion effects. For these indices no 
zero-point corrections have been applied to the SDSS line-strengths as the corrections are not statistically significant, see Section I A4I 
Closed symbols represent indices with SNR > 30 A" 1 , while open symbols have SNR less than this value. The dashed line shows the 
one-to-one correspondence. 



at a ~ 300 km s 1 for the indices that are common between 
our two studies. 



APPENDIX B: TESTING THE MCMC 
APPROACH TO SSP PARAMETER 
ESTIMATION 

To test our MCMC implementation we have chosen to com- 
pare against the Monte Carlo approach of TMBO. We per- 
form a similar test to TMBO, using a model object with 
an age of 10.7 Gyr, [Z/H] = 0.26 and [a/Fe] = 0.25, which 
yields line-strengths of H/3 = 1.59 A, Mgb = 4.73 A and 
(Fe) = 2.84 A. Adopting these line-strengths and the as- 
sociated errors of dU/3 = 0.06 A, dMgb = 0.06 A and 
d (Fe) = 0.05 A we then employ our MCMC approach to re- 
cover the original model object SSP parameters and also 
obtain the errors on the parameters associated with the 
adopted errors on the line-strengths. When evaluating Equa- 
tion [T] we employ uniform base priors which cover the full 
support of each of the SSP parameters and we recover the 
same bivariate parameter distributions as those shown in 
Figure 3 of TMBO. 

We have also estimated the SSP parameters for 123 
early-type galaxies from Table 2 of TMBO and use their 
line-strength data and our MCMC implementation, where 
the MCMC simulation for each object consists of 3 parallel 
chains each of which is 10,000 iterations in size. We use the 
Gelman-Rubin 7?.-statistic to determine when convergence 
to the stationary distribution has occurred and exclude the 
burn-in phase of each chain before pooling the data for the 
three chains and estimating one-dimensional marginalised 
posterior parameters. We show the results of this compari- 
son in Figure iBll 

The results of this comparison indicates that our 
MCMC approach is able to produce reliable parameter esti- 
mates and errors. In Table [Bll our errors are slightly larger 
than those of TMBO but a similar comparison for the data 
in Figure |B"T1 would suggest that our errors are generally in 
good agreement. A more detailed error analysis is prohib- 
ited by the fact that we report errors based on the 15.9 and 
84.2 percentiles of the marginalised posterior distribution, 



Table Bl. Comparing MCMC errors to MC errors from TMBO 



model value TMBO error MCMC error" 



age (Gyr) 


10.7 


±1.48 


+ 1.57 b 

-1.48 


[Z/H] 


0.26 


±0.04 


+0.05 c 
-0.05 


[a/Fe] 


0.25 


±0.02 


+0.05 d 
-0.05 



a errors correspond to the central 68.3% of the distributions 
b base prior, age (Gyr) ~ W(0.31, 18.2) 
c base prior, [Z/H] ~ U{-1.0, 0.97) 
d base prior, [a/Fe] ~ W(-0.25, 0.73) 



and so can obtain asymmetric errors, while TMBO report 
symmetric errors. 



APPENDIX C: LICK INDEX LINE STRENGTH 
DATA 

We provide here the line-strength data for the Lick indices 
used in estimating the SSP parameters in Figures Line 
strengths were measured on coadded spectra calibrated to 
the Lick/IDS system, described in Appendix lAl after select- 
ing individual objects using the strategy described in Sec- 
tion [21 In Tables [CTTIC4I we show line-strengths and their la 
errors. 
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A Age (Gyr) A [Z/H] A [o/Fe] 

Figure Bl. Here we show the residual differences between the MCMC approach of this work and the MC approach of TMBO for 123 
early-type galaxies. For a small number of objects our approach yields slightly older ages than TMBO and consequently lower mctallicitics 
as a result of the age-metallicity degeneracy, hence the tail of the distributions to older ages and lower metallicities. Differences are in 
the sense, this work - literature. 



Table CI. Fully Calibrated Line Strength Measurements for 230 < c v ^ 260 kms 1 mass sample. 



Redshift Interval H/3 ( A) H 7 F ( A) Mgb(A) Fe5270 ( A) Fe5335 ( A) 






02 


< 


z 







04 


1 


728 


± 





046 


-1 


635 


± 





048 


4 


688 


± 





054 


2 


993 


± 





063 


2 


791 


± 





085 





04 


< 


z 







06 


1 


671 


± 





030 


-1 


632 


± 





032 


4 


644 


± 





035 


3 


007 


± 





041 


2 


786 


± 





055 





06 


< 


z 







08 


1 


669 


± 





022 


-1 


600 


± 





023 


4 


669 


± 





027 


2 


956 


± 





030 


2 


743 


± 





039 





08 


< 


z 







10 


1 


655 


± 





023 


-1 


539 


± 





025 


4 


598 


± 





028 


3 


004 


± 





032 


2 


760 


± 





043 





10 


< 


z 







12 


1 


698 


± 





022 


-1 


488 


± 





024 


4 


564 


± 





030 


2 


906 


± 





033 


2 


742 


± 





037 





12 


< 


z 







14 


1 


634 


± 





021 


-1 


484 


± 





022 


4 


549 


± 





029 


2 


917 


± 





027 


2 


752 


± 





031 





14 


< 


z 







16 


1 


660 


± 





025 


-1 


458 


± 





026 


4 


621 


± 





030 


2 


908 


± 





029 


2 


670 


± 





036 





16 


< 


z 







18 


1 


651 


± 





028 


-1 


496 


± 





028 


4 


563 


± 





029 


2 


906 


± 





031 


2 


677 


± 





040 





18 


< 


z 







20 


1 


688 


± 





031 


-1 


371 


± 





030 


4 


559 


± 





030 


2 


843 


± 





033 


2 


623 


± 





040 





20 


< 


z 







22 


1 


654 


± 





034 


-1 


337 


± 





032 


4 


587 


± 





033 


2 


881 


± 





033 


2 


644 


± 





040 





22 


< 


z 







24 


1 


664 


± 





026 


-1 


369 


± 





029 


4 


519 


± 





028 


2 


830 


± 





029 


2 


677 


± 





037 





24 


< 


z 







26 


1 


668 


± 





028 


-1 


268 


± 





034 


4 


530 


± 





030 


2 


843 


± 





033 


2 


535 


± 





039 





26 


< 


z 







28 


1 


698 


± 





032 


-1 


256 


± 





042 


4 


422 


± 





037 


2 


808 


± 





036 


2 


487 


± 





046 





28 


< 


z 







30 


1 


760 


± 





037 


-1 


236 


± 





052 


4 


369 


± 





040 


2 


837 


± 





043 


2 


567 


± 





061 





30 


< 


z 







32 


1 


764 


± 





036 


-1 


126 


± 





058 


4 


285 


± 





043 


2 


793 


± 





052 


2 


554 


± 





062 





32 


< 


z 







34 


1 


772 


± 





042 


-0 


982 


± 





074 


4 


310 


± 





061 


2 


756 


± 





061 


2 


566 


± 





072 





34 


< 


z 







36 


1 


812 


± 





056 


-0 


878 


± 





102 


4 


269 


± 





074 


2 


763 


± 





074 


2 


478 


± 





097 





36 


< 


z 







38 


1 


791 


± 





073 


-0 


701 


± 





129 


4 


005 


± 





101 


2 


768 


± 





119 


2 


456 


± 





185 





38 


< 


z 







40 


1 


885 


± 





121 


-1 


014 


± 





189 


4 


684 


± 





179 


2 


937 


± 





259 


2 


407 


± 





276 
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Table C2. Fully Calibrated Line Strength Measurements for 260 < c v ^ 290kms~ 1 mass sample. 



Redshift Interval H/3 ( A) H7F ( A) Mgb(A) Fe5270 ( A) Fc5335 ( A) 






02 


< 


2 




0. 


04 


1 


.612 


± 


0. 


.052 


-1 


.872 


± 





.054 


4 


.916 


± 





060 


3 


.117 


± 





.071 


2 


.784 


± 





,098 





.04 


< 


2 




0. 


06 


1 


.552 


± 


0. 


031 


-1 


.773 


± 





033 


4 


.934 


± 


0. 


.037 


3 


.054 


± 





.042 


2 


.857 


± 





,059 


0. 


06 


< 


2 




0. 


08 


1 


.659 


± 





.027 


-1 


.603 


± 





028 


4 


.780 


± 


0. 


033 


3 


.014 


± 





.037 


2 


.842 


± 





.050 


0. 


08 


< 


2 




0. 


10 


1 


.519 


± 


0. 


.025 


-1 


.641 


± 





027 


4 


.811 


± 


0. 


032 


3 


.014 


± 





.036 


2 


,771 


± 





,050 


0. 


10 


< 


2 




0. 


12 


1 


.641 


± 





.025 


-1 


.629 


± 


0. 


.027 


4 


.860 


± 


0. 


033 


3 


.004 


± 





.038 


2 


.754 


± 





,044 


0. 


12 


< 


2 




0. 


14 


1 


.582 


± 


0. 


022 


-1 


.578 


± 





.023 


4 


.807 


± 


0. 


031 


2 


.954 


± 





.029 


2 


.797 


± 


0, 


035 


0. 


.14 


< 


2 




0. 


16 


1 


.653 


± 


0. 


024 


-1 


.572 


± 





025 


4 


.785 


± 


0. 


029 


3 


Oil 


± 





.028 


2 


.796 


± 


0, 


036 


0. 


16 


< 


2 




0. 


18 


1 


.646 


± 





.023 


-1 


.556 


± 


0. 


022 


4 


.751 


± 


0. 


.024 


2 


.988 


± 





.025 


2 


.768 


± 





,034 


0. 


18 


< 


2 




0. 


20 


1 


.655 


± 


0. 


.025 


-1 


.432 


± 





023 


4 


.787 


± 


0. 


025 


2 


.980 


± 





.026 


2 


,783 


± 





,033 


0. 


20 


< 


2 




0. 


22 


1 


.725 


± 


0. 


027 


-1 


.396 


± 





.024 


4 


.732 


± 


0. 


027 


2 


.945 


± 





.026 


2 


.834 


± 


0, 


033 


0. 


22 


< 


2 




0. 


24 


1 


.693 


± 


0. 


022 


-1 


.458 


± 





.023 


4 


.696 


± 


0. 


023 


2 


.904 


± 





.023 


2 


,736 


± 


0, 


031 


0. 


2d 


< 


2 




0. 


26 


1 


.680 


± 


0. 


.022 


-1 


.342 


± 





026 


4 


.719 


± 


0. 


.024 


2 


.986 


± 





.026 


2 


.667 


± 





,032 


0. 


26 


< 


2 




0. 


28 


1 


.731 


± 


0. 


026 


-1 


.282 


± 





.033 


4 


.652 


± 


0. 


030 


2 


.916 


± 





.029 


2 


,712 


± 


0, 


038 


0. 


28 


< 


2 




0. 


30 


1 


,755 


± 


0. 


.030 


-1 


.210 


± 





041 


4 


.627 


± 


0. 


032 


2 


.962 


± 





.035 


2 


.763 


± 





,050 


0. 


30 


< 


2 




0. 


32 


1 


.832 


± 


0. 


027 


-1 


.134 


± 





043 


4 


.516 


± 


0. 


033 


2 


.884 


± 





.039 


2 


.655 


± 


0, 


048 


0. 


32 


< 


2 




0. 


34 


1 


.775 


± 





.031 


-1 


.165 


± 





.054 


4 


.484 


± 


0. 


.045 


2 


.900 


± 





,044 


2 


,718 


± 





,055 


0. 


34 


< 


2 







36 


1 


.871 


± 


0. 


.042 


-0 


.982 


± 





077 


4 


.434 


± 


0. 


057 


2 


.864 


± 





,057 


2 


.620 


± 





079 


0. 


36 


< 


2 




0. 


38 


1 


.875 


± 


0. 


.056 


-0 


.999 


± 





099 


4 


.411 


± 


0. 


080 


2 


.926 


± 





,091 


3 


,012 


± 


0, 


146 


0. 


.38 


< 


2 




0. 


40 


1 


.852 


± 


0. 


070 


-0 


.762 


± 





.114 


4 


.461 


± 


0. 


109 


2 


815 


± 





.153 


2 


,715 


± 


0, 


.175 



Table C3. Fully Calibrated Line Strength Measurements for 290 < cr v $5 320kms 1 mass sample. 



Redshift Interval H/3 ( A) H-yF ( A) Mgb(A) Fe5270 ( A) Fe5335 ( A) 






02 


< 


2 







04 


1 


.534 


± 





,073 


-1 


.675 


± 





.074 


5, 


.055 


± 





.087 


3 


,115 


± 





,103 


2 


,771 


± 





,150 





04 


< 


2 




0. 


06 


1 


,481 


± 





082 


-1 


.918 


± 





.088 


5. 


160 


± 





.106 


3 


.143 


± 





.121 


3 


,107 


± 





,171 


0. 


06 


< 


2 




0. 


08 


1 


.524 


± 


0, 


044 


-1 


.861 


± 





,048 


5. 


197 


± 





.057 


3 


.040 


± 





.065 


2 


,971 


± 





,090 


0. 


08 


< 


2 




0. 


10 


1 


.566 


± 


0, 


,042 


-1 


.781 


± 





,047 


5. 


008 


± 





,056 


3 


.083 


± 





,063 


2 


.961 


± 





.090 


0. 


10 


< 


2 




0. 


12 


1 


.612 


± 





,036 


-1 


.705 


± 





040 


5. 


058 


± 





.051 


3 


.083 


± 





,057 


2 


.899 


± 


0. 


,070 


0. 


12 


< 


2 




0. 


14 


1 


.586 


± 





,030 


-1 


.684 


± 


0. 


,033 


4, 


963 


± 





.045 


3 


039 


± 





,041 


2 


.871 


± 


0. 


.051 


0. 


14 


< 


2 




0. 


16 


1 


.567 


± 





029 


-1 


.565 


± 





.030 


5. 


016 


± 





.037 


3 


Oil 


± 





,035 


2 


.808 


± 


0, 


.048 


0. 


16 


< 


2 




0. 


18 


1 


.630 


± 


0, 


,031 


-1 


.623 


± 





031 


4 


.935 


± 





.035 


3 


,040 


± 





,036 


2 


.840 


± 


0. 


051 


0. 


18 


< 


2 




0. 


.20 


1 


.626 


± 





,029 


-1 


.548 


± 





027 


4, 


997 


± 





.030 


3 


110 


± 





,032 


2 


.843 


± 





,041 


0. 


20 


< 


2 




0. 


22 


1 


.622 


± 


0, 


.029 


-1 


.467 


± 





.026 


4. 


,943 


± 





.030 


3 


.075 


± 





,029 


2 


.851 


± 


0, 


.038 


0. 


22 


< 


2 




0, 


24 


1 


.640 


± 


0, 


025 


-1 


.505 


± 





026 


4. 


.896 


± 





.027 


2 


,990 


± 





,027 


2 


.746 


± 


0, 


038 


0, 


24 


< 


2 




0. 


.26 


1 


.645 


± 





,024 


-1 


.423 


± 





,028 


4. 


845 


± 





.027 


3 


.060 


± 





,029 


2 


.799 


± 





,037 


0, 


26 


< 


2 




0. 


28 


1 


.728 


± 





,026 


-1 


.376 


± 





033 


4, 


913 


± 





.031 


3 


,022 


± 





,029 


2 


.799 


± 





,040 


0, 


28 


< 


2 




0. 


30 


1 


.755 


± 


0, 


027 


-1 


.291 


± 





,037 


4 


735 


± 





.030 


3 


015 


± 





,031 


2 


.803 


± 


0, 


048 


0, 


30 


< 


2 




0. 


32 


1 


.759 


± 


0, 


029 


-1 


.269 


± 





.045 


4. 


733 


± 





.035 


3 


.037 


± 





.042 


2 


.796 


± 


0, 


.054 


0. 


32 


< 


2 




0, 


34 


1 


,781 


± 





,031 


-1 


.087 


± 





052 


4. 


635 


± 





.045 


2 


.997 


± 





,044 


2 


.629 


± 





,058 


0. 


34 


< 


2 




0. 


.36 


1 


.805 


± 





,044 


-1 


.159 


± 





,080 


4 


546 


± 





.062 


2 


837 


± 





.060 


2 


.847 


± 


0, 


.088 


0. 


36 


< 


2 




0. 


38 


1 


.779 


± 


0, 


047 


-1 


.170 


± 





.082 


4. 


457 


± 





.069 


3 


,089 


± 





,077 


2 


.533 


± 


0, 


131 


0. 


,38 


< 


2 




0. 


40 


1 


.805 


± 





.062 


-1 


.155 


± 





100 


4. 


577 


± 





.097 


2 


.869 


± 





,136 


2 


.557 


± 


0, 


166 
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Table C4. Fully Calibrated Line Strength Measurements for 320 < c v ^ 350kms _1 mass sample. 



Redshift Interval H/3 ( A) H7F ( A) Mgb(A) Fe5270 ( A) Fc5335 ( A) 






02 


< 


2 




0. 


04 


1 


.406 


± 





139 


-1 


.820 


± 





.134 


5. 


.270 


± 


0. 


170 


3 


.314 


± 


0. 


.194 


3 


.061 


± 





.287 





.04 


< 


2 




0. 


06 


1 


.484 


± 


0. 


086 


-1 


.893 


± 





087 


5. 


132 


± 





109 


3 


.197 


± 


0. 


123 


3 


.000 


± 





.185 


0. 


06 


< 


2 




0. 


08 


1 


,481 


± 





.109 


-1 


.717 


± 





.123 


5. 


405 


± 


0. 


148 


3 


.316 


± 


0. 


161 


3 


.053 


± 





.239 


0. 


08 


< 


2 




0. 


10 


1 


.518 


± 





.067 


-1 


.895 


± 





076 


5. 


359 


± 


0. 


094 


3 


.064 


± 





105 


3 


.038 


± 





.151 


0. 


10 


< 


2 




0. 


12 


1 


.662 


± 





.075 


-1 


.794 


± 





.087 


5. 


147 


± 


0. 


113 


3 


107 


± 


0. 


124 


3 


.130 


± 





.157 


0. 


12 


< 


2 




0. 


14 


1 


.504 


± 





059 


-1 


.658 


± 





.066 


5. 


122 


± 


0. 


094 


3 


133 


± 


0. 


084 


2 


.766 


± 


0. 


.110 


0. 


.14 


< 


2 




0. 


16 


1 


.658 


± 


0. 


053 


-1 


.673 


± 





.056 


5. 


.229 


± 


0. 


070 


3 


126 


± 


0. 


067 


3 


.019 


± 


0. 


095 


0. 


16 


< 


2 




0. 


18 


1 


.651 


± 





.047 


-1 


.713 


± 





.049 


5. 


153 


± 


0. 


056 


3 


100 


± 


0. 


058 


3 


.067 


± 





.084 


0. 


18 


< 


2 




0. 


20 


1 


.649 


± 





.043 


-1 


.591 


± 





041 


5. 


181 


± 


0. 


048 


3 


.131 


± 


0. 


050 


3 


.051 


± 





.068 


0. 


20 


< 


2 




0. 


22 


1 


.603 


± 


0. 


043 


-1 


.544 
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APPENDIX D: DERIVED SSP PARAMETER 
DATA 

We provide here the SSP parameters determined for each of 
the mass samples shown in Figure [6] Parameter estimates 
are determined from the posterior probability distribution 
described in Section [3] and Section [3] Briefly, each SSP 
parameter is estimated by determining the 16, 50 and 84 
percentiles of posterior distribution after marginalising over 
the remaining two parameters. For Tables IDUD4I the pos- 
terior distribution has been obtained using uniform priors 
that span the full support of each of the KMT05 parame- 
ters, i.e. age (Gyr) ~ W(0.31, 18.2), [Z/H] ~ W(-1.0, 0.97), 
[a/Fe] ~ W(-Q.25, 0.73). For Table EUthe posterior distri- 
bution has been re-weighted, as described in Section 21 to 
account for Gaussian priors placed on [Z/H] and [a/Fe] , in 
order to better constrain the age-redshift relationship. 
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Table Dl. SSP parameter data for 230 < a v ^ 260 km s 1 mass Table D2. SSP parameter data for 260 < a v ^ 290 km s 1 mass 

sample using H7F, Mgb , Fe5270, Fe5335. sample using H7F, Mgb , Fe5270, Fe5335. 



Redshift Interval Age (Gyr) [Z/H] [o/Fe] Redshift Interval Age (Gyr) [Z/H] [a/Fe] 
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Table D3. SSP parameter data for 290 < u v t< 320 km s" 1 mass 
sample using H7F, Mgb , Fe5270, Fe5335. 



Redshift Interval Age (Gyr) [Z/H] [a/Fe] 
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Table D4. SSP parameter data for 320 < a v ^ 350 km s 1 mass 
sample using H7F, Mgb , Fe5270, Fe5335. 



Redshift Interval Age (Gyr) [Z/H] [a/Fe] 
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Table D5. SSP ages determined for all mass samples using H7F, Mgb , Fe5270, Fe5335 
and using Gaussian priors on [Z/H] and [cc/Fe] . 



Redshift Age (Gyr) Age (Gyr) Age (Gyr) Age (Gyr) 

Interval 230< a v <260 260< a v ^290 290< a v ^320 320< a v SC350 
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